Polarization field in a single-valley strongly-interacting 2D electron system 
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The magnetic field of complete spin polarization is calculated in a disorderless single-valley 
strongly-interacting 2D electron system. In the metallic region above the Wigner-Mott transition, 
non-equilibrium spin states are predicted, which should give rise to hysteresis in the magnetization. 
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A transition from the high conductivities er 3> e 2 //i 
to the low conductivities a <C e 2 //i that is observed 
with decreasing electron density in disordered 2D elec- 
tron systems was first interpreted as the Anderson metal- 
insulator transition [l[ . The interpretation changed after 
the appearance of the weak localization theory [2[ and 
scaling hypothesis Q. According to the latter, there is 
no metallic state in disordered infinite 2D electron sys- 
tems at zero temperature. As a result, the transition 
observed in experiment was referred to as "apparent" 
metal-insulator transition. 

The next interpretation change was caused by the find- 
ing that the temperature dependence of the conductivity 
in low-disordered 2D electron systems changes sign at a 
critical density attributed to the metal-insulator transi- 
tion point The observed critical density corresponds 
to strong electron-electron interactions, whereas nonin- 
teracting or weakly-interacting 2D systems are consid- 
ered in the scaling concept. Moreover, the spin suscepti- 
bility and effective mass in the least-disordered 2D elec- 
tron systems were found to increase at low densities with 
a tendency to diverge at some electron density that is 
weakly dependent on disorder [f| Q . This indicates that 
the metal-insulator transition observed in a number of 
2D electron systems is driven by interactions [7j, while 
the disorder is of minor importance. Additional confir- 
mation for this statement is provided by the observation 
of the critical increase of the effective mass with increas- 
ing interactions in the 2D fermion system composed of 
He 3 atoms @]. 

There has been published a good deal of experimen- 
tal and theoretical work on the metal-insulator transi- 
tion in two dimensions and related phenomena (see, e.g., 
experimental reviews [^4l2]| and theoretical publications 
[T3l - tl7| ). The key question is whether the metal- insulator 
transition observed in low-disordered 2D electron systems 
is a transition to the Wigner crystal. The mass diver- 
gence when the crystallization point is approached on 
the metallic side was already demonstrated in the paper 
[l8| using Gutzwiller's theory. Recently, many properties 
of strongly-interacting 2D and 3D systems have been ex- 
plained within the concept of the Wigner-Mott transition 
using dynamical mean-field theory |17| . 

In this paper, we study the magnetic field of full spin 
polarization B c as a function of electron density n s for 
the Wigner-Mott transition. The magnetic field is as- 
sumed to be parallel to the 2D electron system and act 



on electron spins only. We use a lattice model [T8j by 
adapting the results of the papers [H, Ho| to the problem 
to be solved. The real 2D electron system is replaced 
by lattice sites with density n s . If for the lattice in- 
troduced there is an energy gap between the lowest and 
next bands, filling every site with one electron results in a 
Wigner crystal model. However, in the ground state, two 
electrons with anti-parallel spins can be located, despite 
repulsion, on each site with some probability depend- 
ing on n s . Such sites determine the portion of mobile 
electrons and, thus, the transport properties of the elec- 
tron system. For the case of noninteracting electrons, 
the model used corresponds to the half-filled lowest en- 
ergy band. The noninteracting electrons in zero mag- 
netic field occupy all states up to the Fermi momentum 
kp, each state being filled with spin- up and spin-down 
electrons. In the lattice model, the probability to find on 
each site two electrons is equal to 1/4, which is the same 
as the probability to find on each site no electrons. 

The momentum distribution function for the system 
with electron interactions in the metallic state is shown in 
Fig. [lja) for zero magnetic field. All states are occupied 
up to the Brillouin zone boundary 1zq. A jump of the dis- 
tribution function with height Z < 1 occurs at kp . With 
decreasing electron density both kg and kp decrease, the 
ratio kp/ko being constant. As the critical density n c is 
approached, the jump Z — > and the distribution func- 
tion f ud — >• 0.5. The shaded region that spans up to 
ko in the figure corresponds to the Wigner crystal. The 
many-electron wave function [n| describes the superpo- 
sition of the crystal state and the Fermi quasi-particles 
with an effective mass m* that is inversely proportional 
to the distribution function jump at the Fermi momen- 
tum kp: m* /m — 1/Z, where m is the band mass. 

The distribution function for electrons with different 
spins in magnetic fields is displayed in Fig. [TJb, c). Both 
the quasi-particle weight Z and the Fermi momenta of 
quasi-particles with different spins are influenced by mag- 
netic field. Evidently, the increase of magnetic field leads 
to k F — > ko and k F — > 0. In what follows we find that 
with increasing magnetic field the quasi-particle weight Z 
increases (decreases) for n s 3> n c (n s > n c ) until at some 
degree of spin polarization, the electron system reaches 
by jump full spin polarization in which case every lattice 
site is filled with one electron. The presence of the jump 
of the degree of spin polarization gives rise to hysteresis 
with decreasing magnetic field. 
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FIG. 1: Probability to occupy a state with momentum k as 
a function of momentum in a certain direction: (a) for spin- 
unpolarized system; (b) for spin-up electrons; (c) for spin- 
down electrons. The shaded region corresponds to the Wigner 
crystal. 



In accordance with the results [ljj, [2l|, the ground 
state-averaged Hamiltonian in a parallel magnetic field 
B is written 



<H>- 



4Dn 



where v is the portion of doubly occupied sites (0 < v < 
1/2), e is the electron charge, K is the dielectric con- 
stant, Do is the density of states of the spin-polarized 
noninteracting electrons. The degree of spin polarization 
P = (n u — rid)/n s is determined by the difference of the 
densities of spin-up and spin-down electrons, and the co- 
efficient a is determined by the electron wave function on 
the lattice site. The kinetic energy is counted from the 
band-averaged energy for noninteracting electrons. The 
relation between Z and v taking account of correlations 
is as follows 
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For p = 0, Eqs. (Q] [2]) reduce to the known ones in zero 
magnetic field |18j |. 

It is necessary to find a minimum of the expression ([1]) 
over v, regarding the relation ([2]): 



d < H(v,p,n s ) > 
dv 

which yields a dependence v{jp,n s ) 

(1 - 2v)(l - Av) - p 2 
[(l-2i/) 2 -p 2 ] 1 /2 



0. 



(3) 
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where n c = (eve 2 _D /2k) 2 . 

Note that for p = 0, the ground state energy is a min- 
imum if 



1/2' 
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FIG. 2: (a) Portion of doubly occupied sites and (b) the quasi- 
particle weight versus spin polarization at different values of 
the interaction parameter. 



At n s — > n c , both v and Z tend to zero, i.e., every lattice 
site is filled with one electron, and the quasi-particle mass 
diverges 
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In the opposite limiting case of n s — > oo, the values v and 
Z are v = 1/4 and Z = 1, and the quasi-particle mass is 
equal to to* = to. 

The dependence v{p) is shown in Fig. HJa) for different 
values of the interaction parameter n c /n s . If n s > 4ri c , 
the value v zeroes only at p = 1, whereas at n s < 4ra c 
the equality v — is the case over a whole range of p. 
We note that this fact is important for determining the 
dependence B c (n s ), and the density n s = 4n c manifests 
itself as "apparent" critical density. The corresponding 
dependence Z(p) is represented in Fig.^b). The value Z 
increases (decreases) with spin polarization at n s > 4n c 
(n s < 4n c ), attaining the limiting value Z = 1 [Z = 0). 
Interestingly, the spin states near p — 1 turn out to be 
non-equilibrium states. If the dependence < H(p) > is 
two minima separated by a maximum, the electron sys- 
tem at T = gets to the deeper minimum in a magnetic 
field where the maximum disappears, which causes hys- 
teresis. 

We shall start by considering the high- density region 
n s > 16n c . In the spirit of the paper [22j, one needs to 
find a minimum of the Hamiltonian (JTJ [2]) over p: 



d < H(i/,p,n s ) > 



dp 



= o, 



(4) which yields 
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By solving Eqs. ((HIS]) with p — > 1, we get the polarization 
field 



1 - 2 



1/2- 



(9) 



3 



-0.01 



Q 

A 

a 

V 



-0.02 




-0.03 



FIG. 3: The ground state energy as a function of spin polar- 
ization at fixed n c /n s = 0.6 for different magnetic fields. The 
minimum is marked by arrows. 



The dependence obtained is the same as the one sug- 
gested earlier [5[ and is similar to B c (n s ) observed in 
experiment. Obviously, Eq. (j9]) cannot describe the equi- 
librium value B c down to n s = 4n c , because the Wigner- 
Mott transition in zero magnetic field occurs at the lower 
density n s = n c . 

We solve the problem numerically at electron densities 
near 4n c . There, the system behavior turns out to bc 
hysteretic. Using the dependence v(p) of Fig. [3J we cal- 
culate the dependence < H(p) > for a fixed value of the 
interaction parameter n c /n s at different magnetic fields, 
as shown in Fig. [3] for n c /n s — 0.6. In weak magnetic 
fields, the energy has a minimum corresponding to the 
realized degree of spin polarization. As the magnetic 
field is increased, the minimum turns into an inflection 
and disappears entirely. The inflection determines the 
upper polarization field at which the system state jumps 
to p = 1, v = 0. The degree of spin polarization first 
increases linearly with magnetic field, then turns up, and 
finally, the electron system reaches complete spin polar- 
ization by jump, as shown in Fig. |4j With decreasing 
magnetic field there arises hysteresis. The electron sys- 
tem stays in the p — 1 state until the maximum at p < 1 
in the dependence < H(p) > (Fig. [3]) disappears at the 
lower polarization field at which the system state jumps 
out of p = 1 to close the hysteresis loop (Fig. @| . It is 
clear that the lower polarization field is given by Eq. © 
at n s > 4n c and is equal to zero at n s < 4n c . 

The dependence of the polarization field on electron 
density with increasing and decreasing magnetic field is 
displayed in Fig. [5] along with the effective mass as a 
function of electron density in zero magnetic field. Based 
on linear extrapolation of the high-density values, the 
field B c tends to zero at n s 4n c , which is well above 
the density where the effective mass in zero magnetic field 
diverges. Note that very low temperatures are needed to 
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FIG. 4: Hysteresis of the degree of spin polarization with 
changing magnetic field at fixed (a) n c /n s = 0.6 and (b) 
n c /n„ — 0.2. 



observe the hysteresis in B c because the maximum in 
Fig. |31 is smeared due to thermal fluctuations already at 
T - 0.1 K. 

In accordance with the solution, one might conclude 
that the magnetic field promotes crystallization at n s < 
An Cl whereas at higher densities, the system state in 
magnetic fields is the electron liquid. As a matter of 
fact, the complete filling of the first Brillouin zone in 
the model used means crystallization only if the energy 
gap is present for the model lattice. In turn, the en- 
ergy gap is determined by electron-electron interactions 
on different lattice sites so that in the frames of the model 
[l8l - l20l ] , one cannot say whether or not the energy gap is 
present. Nevertheless, one can expect that the predic- 
tion is qualitatively correct: the gap in magnetic fields 
should survive at electron densities above n c and disap- 
pear at yet higher densities. Within the hysteresis un- 




FIG. 5: Polarization field as a function of electron density 
with increasing (solid black line) and decreasing (dashed red 
line) magnetic field. Also shown by the dash-dotted blue line 
is the effective mass in B = versus electron density. 
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certainty, the low-density part of B c (n s ) separates the 
correlated electron liquid and crystal and corresponds to 
the metal-insulator transition. The last fact is confirmed 
by comparison of the experimental results of Ref. Q and 
Refs. dlUl. 

Another prediction of our calculations is that the ef- 
fective mass in a single- valley electron system should de- 
crease with magnetic field and reach the noninteracting 
electron mass at the onset of complete spin polarization if 
the system stays in the metallic state. The same behavior 
of the mass should be the case for full isospin polarization 
in a two-valley spin-polarized electron system. Accord- 
ing to the statements of the papers [IH Hf|, the effect 
has been observed recently in the 2D electron system in 
AlAs quantum wells. 

The magnetic field of full spin polarization as a func- 
tion of electron density in the single- valley 2D electron 
system of GaAs/AlGaAs heterostructures was studied in 
Ref. [27|. Because of large thickness of the 2D system, 
the dependence B c (n s ) is strongly distorted by orbital 
effects. Nevertheless, its slope for low B c is in agree- 
ment with the slope expected from Eq. 0, and B c (n s ) 
extrapolates to zero at a finite density [lfj . 

The most complete and detailed experimental infor- 
mation has been obtained for two- valley 2D electron sys- 
tems. However, our calculations are not relevant for such 
systems, and comparison with those experiments is not 
justified. 

It is interesting to compare our results with alterna- 



tive calculations. The behavior of the effective mass near 
the critical density of Eq. © has been reproduced in 
Ref. [13] • There, the dependence B c (n s ) has also been 
obtained which, unlike our results, does not reveal hys- 
teresis. 

It is worth noting that the polarization field B c versus 
n s for a two- valley 2D electron system has been calcu- 
lated using quantum Monte-Carlo simulations (22[. In 
the clean limit (kpl 3> I, where I is the mean free path), 
the divergence of the mass and the critical B c (n s ) are 
not found. However, it is strange that the degree of spin 
polarization in the paper [22| is proportional to magnetic 
field, regardless of disorder and interaction strength. 

In summary, we have calculated the dependence of the 
magnetic field of complete spin polarization on electron 
density for the Wigner-Mott transition in a single-valley 
2D electron system. The following predictions of the 
model used have not yet been confirmed in experiment. 
In the metallic region above the Wigner-Mott transition, 
non-equilibrium spin states are expected, which should 
lead to a hysteretic behavior of the magnetization. Lin- 
ear extrapolation to zero of the high-density values of the 
polarization field yields a critical density that exceeds the 
Wigner-Mott transition point by the factor of about four. 
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Ministry of Sciences. 
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